Charge fluctuations and the tunneling spectra of non-magnetic metallic nanoparticles 
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We present microscopic transport calculations of the tunneling spectra of non-magnetic metal 
nanoparticles. We show that charge fluctuations give rise to tunneling resonances of a new type. 
Positive and negative fluctuations have differing kinetics and thus account for previously unexplained 
spectral features that are found experimentally under only forward or only reverse applied bias. The 
observed clustering of tunneling resonances of Al nanoparticles arises naturally from our theory. 

PACS numbers: 73.22.-f, 73.22.Dj 
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In metal particles with dimensions on the nanometer 
scale the electron de Broglie wavelength is comparable to 
the size of the particle. Thus these nanoparticles exhibit 
discrete electronic spectra that can be observed directly 
in tunneling measurements. |l[] Such experiments have 
recently been carried out for Al, [|| Au, Ag, Cu jj) and 
Co |Q nanoparticles coated with aluminum oxide that 
forms the tunnel barrier. They have attracted consid- 
erable attention since they can in principle provide de- 
tailed microscopic information relevant to many impor- 
tant but poorly understood aspects of nanoscale metal 
physics that range from the effects of disorder and surface 
chemistry to nanoscale ferromagnetism and superconduc- 
tivity. [I]] Some of the observed effects have been mod- 
elled phenomenologically with considerable success. |^-||] 
However even in the case of Al nanoparticles (the sim- 
plest and most studied of these systems) the present un- 
derstanding of the results of the experiments is far from 
satisfactory. For example, some of the tunneling reso- 
nances that are seen experimentally can be matched with 
similar features that are observed when the bias applied 
to the nanoparticle is reversed. Thus they can reasonably 
be attributed to tunneling through particular electronic 
states of the nanoparticle. || However, other observed 
tunneling resonances have no identifiable counterparts 
under reverse bias and their physical origin has remained 
a mystery. This suggests that some important aspects 
of the physics of electron transport through the metal 
nanoparticles have not been recognized to date. In this 
Communication we identify a plausible candidate: We 
demonstrate theoretically that charge fluctuations that 
occur whenever a current flows through the nanoparti- 
cle should result in transport resonances of a new type 
that begin already in the first step of the Coulomb stair- 
case. We predict that for most samples these new reso- 
nances (unlike other tunneling features) should be much 
stronger for one direction of the applied bias than for 
the other. Thus charge fluctuations account naturally for 
the presence in the experimental tunneling spectra of the 
previously unexplained resonances described above. We 
argue that in typical samples this new mechanism should 
account for a substantial fraction of all of the observed 
tunneling features. 



We illustrate our predictions with numerical calcu- 
lations for Al nanoparticles whose electronic structure 
is described by a microscopic tight-binding model [Tc| ] 
that incorporates the geometry of the particle and ac- 
counts for the presence of disorder as well as the de- 
tailed chemistry of the metal-oxide interface. Salient fea- 
tures of the electronic structure are depicted in Fig. |l| 
which shows the calculated energy eigenvalues Ei near 
the Fermi level Ep of a disc-shaped nanoparticle of vol- 
ume V = 16.9nm 3 . The amplitudes of the electron eigen- 
functions ipi on the top (T) and bottom (B) surfaces of 
the nanoparticle are also displayed for selected electron 
eigenstates. The diameter of each circle represents the 
magnitude of ipi at a given atomic site (indicated by the 
central dot). Due to the presence of surface disorder, 
the relief of the amplitude is quite complex and depends 
strongly on the electronic state that is considered and on 
which surface is shown. 

The wavefunction landscapes at these surfaces enter 
the electron tunneling efficiencies 7^ between level i 
of the nanoparticle and contact A (—T,B) as follows: 
Let MA be the tunneling matrix element between the 

electronic state \&e(r) = (f*|^e) °f the contact and 
the state fi{r) = (r\ifi) of the nanoparticle. Since 
MA is a transfer matrix element we adopt the expres- 
sion MA = 1Z J n ^dS ty^ifi widely used in the quan- 
tum chemistry literature to calculate transfer matrix el- 
ements between different molecular states. 11 1Z is an 
energy scale factor and 0> is the T or B surface of the 
nanoparticle. The nanoparticle electron wavefunction is 
fiif) = J2r . a< j.i4>a(r — Rj) where the coefficients a" fi 
arise from diagonalization of a tight-binding Hamiltonian 
constructed using a s, p, d orbital basis (<f> a ). |10| Ob- 
taining a realistic fj(r) in the contact /oxide region is 
difficult. However, we do not need the whole wavefunc- 
tion but only the values that ^f^(r) takes at the surface 
of the nanoparticle. From a tight-binding point of view 
with only s orbitals in the basis, the wavefunction at 
the surface can be written ^^(r) = J^p ^j,e^s(^ — ^i) 
with Tj the positions of the atomic sites at VL\. Due to 
the disordered nature of the metal-oxide interface and 
time reversal symmetry, we represent the coefficients bj e 
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by those arising from diagonalization of a matrix M. 
of order M within the Gaussian Orthogonal Ensemble 
weighted by the WKB extinction coefficient 



1- 



Af is the number of sites at surface k = y2mVb/h 
and Vb ~ 1.2eV |p^| ; d°^° is the thickness of the oxide 
barrier. We average over an ensemble of matrices A4 to 
model the fact that different lead states |^ e ) couple to a 
given nanoparticle state \tpi). Thus we find 
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v\{Eh) is the density of states at the highest occupied 
level in the lead A and (• • •} denotes ensemble averaging. 

In Eq. [I] 7* depends on the oxide thickness be- 
tween the A-lead and the nanoparticle. Within a simple 
parallel-plate capacitor model this thickness also deter- 
mines the lead-dot capacitances Ct and Cb for leads T 
and B; for simplicity we assume Crd%? = Csd^. Fig. || 
shows the calculated 7* , in units of T = 2-k1Z 2 v\{E h ) /K 
for the T and B surfaces, for the first few energy levels 
around the Fermi level for different capacitance ratios 
(Cb/Ct — 1, 1.25). We take as a reference dg" = 5 A and 
k = 0.56A _1 EJ. When the thicknesses (capacitances) 



are equal, the tunneling efficiencies, although asymmet- 
ric, are of the same order of magnitude. As soon as the 
thickness (capacitance) ratio is changed slightly so that 
d°B < dp* , the asymmetry becomes stronger and falls 
rapidly by almost an order of magnitude relative to 7^ . 

We calculate the electric current I through the 
nanoparticle and the differential conductance dl/dV as- 
suming for simplicity that when an electron enters or 
leaves the nanoparticle the latter relaxes quickly to its 
electronic ground state. The complementary limit of 
slow relaxation has been studied previously (|J^], omit- 
ting however the new effects of charge fluctuations that 
we introduce here. We consider forward bias (FB) and 
reverse bias (RB) voltages V at which the electron popu- 
lation N of the nanoparticle changes by no more than ±1 
from its neutral value Nq. Thus the nanoparticle is neu- 
tral (n — N — Nq = 0) or negatively (n = 1) or positively 
(n = —1) charged. We treat the electrostatic charging en- 
ergy of the nanoparticle and leads according to standard 
Coulomb blockade theory [jl5| . We assume that at V = 
the Fermi levels of the leads align with the highest oc- 
cupied level of the neutral particle. The electrochemical 
potential of the T(B) lead is /i T(S) = E F ± eVC s(T) /C s 
with Cs = Ct + Cb [f[6|- The master equation for the 
population of the nanoparticle is then 
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where the Ei are the one-electron energy levels of the 
neutral nanoparticle, U = e 2 /2Cs is the single-electron 
charging energy Jig] , fix) is the Fermi distribution of 
the leads at temperature Td and 9{x) is the Heaviside 
function describing the occupation of the levels in the 
nanoparticle. E F is the highest occupied level with spin 

er(=f, J.). Ep^" and E F +(7 denote the highest occupied 
level of the nanoparticle with spin a after removal(— ) or 
addition(+) of an electron with spin a'. In Eq. we 
assume FB with electron flow from T to 5, therefore 
fi T > n B . In RB the indices T and B are interchanged. 

We now discuss the physical meaning of the terms in 
Eq. ||, keeping for the sake of clarity to the case Td = 0. 
The sum over i in Eq. || describes processes in which an 
electron moves from lead T to the single-electron state 
\ipi) of an initially neutral (n = 0) nanoparticle. This is 
allowed energetically if Ei + U < Ep + eCsV/C^ 1 i.e., 
the single-electron energy Ei of the level in the dot to 
which the transition is made plus the charging energy 
U must be lower than fj, T , the highest occupied level in 
the electrode. Also the single-particle state \(fi) of the 



nanoparticle that accepts the electron, must be initially 
unoccupied. f{x) and 9(x) with the arguments in Eq. |2| 
account for these constraints. Analogous reasoning leads 
to the other terms in Eq. || The sum over j describes 
tunneling from the neutral particle to electrode B which 
is allowed energetically if Ej~U > /i S . If the nanoparti- 
cle is initially charged different energetic restrictions ap- 
ply: If n = —1 initially, tunneling from electrode T[B) to 
the nanoparticle state \ipk) is allowed if E k — U < /j, t ^ b >. 
If n = 1 initially, an electron may tunnel from state \tpi) 
of the nanoparticle to electrode B(T) if Ei + U > fx B ^ T \ 

The rate equations for the probability of the nanopar- 
ticle being in a given charge state (n = — 1,0, +1) are: 
8 t Pi = S1P0 - S4P1, and dtP-i = S 2 P - S3P-1, with 
Po + Pi + P_i = 1. In forward bias the total rates are 
S 1 = S^ d , S 2 = S^ B , S 3 = Sj- d + S B ^ d , and 
£>4 = S d ^ B + Sf^ T . We solve these rate equations for 
the steady state occupation probabilities P**. Then the 
current I passing through the nanoparticle is 

X FB = e(S^ d P? + Si- d P!_\ - S d ^ T Pf). (3) 
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Figure |3j shows our calculated differential conductance 
spectra for Al nanoparticles within the first step of the 
Coulomb staircase. The resonances labelled iX (Xi) are 
due to tunneling from (to) level \Ep + i) of the neutral 
nanoparticle to (from) contact X — T,B. Here i = —i. 
We refer to these resonances as neutral peaks (NP). The 
other peaks labelled Q + {Q) are due to the n = 1 (—1) 
charge fluctuations on the nanoparticle that are induced 
by the current. They occur when a level of the charged 
nanoparticle (renormalized by the charging energy U) be- 
comes available for transport as the electrochemical po- 
tentials of the contacts are swept. These resonances will 
be referred to as charge fluctuation peaks (CFP). Since on 
the first step of the Coulomb staircase the nanoparticle 
is neutral most of the time, most NP's are much stronger 
than most CFP's. However as will be explained below, 
some CFP's can be strongly enhanced by kinetic bottle- 
necks and thus become very prominent spectral features. 

Another striking aspect of Fig. |^ is the strong ten- 
dency of the resonances to cluster that is also found ex- 
perimentally Q. To account for this clustering, previous 
theories employed the phenomenological assumption 
that excited nanoparticle states are very long lived. Here 
we do not make this assumption; the clustering in Fig. || 
follows directly from our microscopic theory |l0| of the 
electronic structure of the nanoparticle. Also in common 
with the experimental data (but unlike the results of 
previous theories || ) the amplitudes of the peaks in Fig. 
H show no systematic decrease with increasing bias V. 

We now discuss the tunneling spectra in Fig. || in more 
detail. For FB we assume n T > fi B so the electrons flow 
from T to B. Note that the highest occupied level of the 
neutral nanoparticle is i = and it is doubly occupied. 

In Fig. ||a Cb — Ct so the applied bias is divided 
equally between the two contacts. In FB the first reso- 
nance is OB (tunneling from the i — neutral nanopar- 
ticle level to lead B) and occurs at V = 50mV= 2U/e. 
This is followed at higher bias by other NP's interspersed 
with CFP's Q + and Q. Since Cb — Ct the spectrum un- 
der RB is similar with the roles of contacts T and B 
interchanged; each NP and CFP under FB has a coun- 
terpart at the same voltage under RB. The amplitudes of 
corresponding peaks under FB and RB are not identical 
because the calculated jf' B are unequal; see Fig. ^a. 

Experimentally, the B and T tunnel barriers have dif- 
ferent thicknesses and hence unequal capacitances and 
very different tunneling efficiencies; see Fig. |^b. Fig. 
|3|b shows such a case with Cb — 1.25Ct which implies 
bias-induced changes in (x T 1.25 times larger than those 
in [i B . Thus corresponding NP's now occur at FB and 
RB values that differ by the capacitance ratio. For ex- 
ample, OT is at 45mV in RB while its FB partner OB is 
at 45mVxl.25 =56mV. Despite the large asymmetry be- 
tween and 7 4 T in Fig. |^b, the amplitudes of the NP's 
are quite similar under FB and RB in Fig. ||b because 
in each case the current must pass through both the T 
and B tunnel barriers. However for asymmetric tunnel 



barriers the amplitudes of the CFP's are very different 
under FB and RB: In Fig. |b the Q CFP's are much 
stronger than the Q + CFP's under FB while the reverse 
is true under RB. This is because asymmetric tunnel bar- 
riers affect the kinetics of positive and negative charge 
fluctuations differently: A Q + charge fluctuation implies 
an electron surplus on the nanoparticle. Under FB this 
fluctuation dissipates easily to the drain (the B contact) 
since is large, thus the total current and dl/dV are 
affected little by the introduction of an additional decay 
channel and the Q + CFP's are weak. Conversely, a Q 
fluctuation (an electron deficit on the nanoparticle) does 
not dissipate as readily under FB since this requires tun- 
neling from the source (T) contact through the weakly 
transmitting (T) barrier. Thus here the introduction of 
a new conducting channel has a larger effect on / and 
dl/dV and Q CFP's are strong in FB. Under RB the 
transport bottleneck is reversed and Q CFP's are weak 
while Q + CFP's are strong. Thus the new charge fluc- 
tuation resonances that we have introduced here readily 
account for the previously puzzling experimental tunnel- 
ing features that have no identifiable partner when the 
bias is reversed. 

In conclusion, we have presented the first calculations 
of electron transport through metal nanoparticles that 
are based on a microscopic theory of their electronic 
structure. We have shown theoretically that charge fluc- 
tuations give rise to tunneling resonances of a new type 
that account for the behavior of previously unexplained 
spectral features that are observed experimentally. 
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FIG. 1. Energy levels near Ef (upper panel). Amplitude of ifi at the T and B surface corresponding to the second level above 
(right) and below (left) the Fermi level (center). The second and third levels above the Fermi energy are nearly degenerate on 
this scale. 5 Al = (4E Al /SAf^V' 1 = 5meV. Ef is the bulk Fermi energy and N Al the electron density of Al. 



FIG. 2. Tunneling efficiencies for energy levels \Ef +i) (i = —i) for T and B leads. 

FIG. 3. dl/dV in forward (FB) and reverse (RB) bias on a logarithmic scale; consecutive ticks indicate a decade. U = 25 
meV, T d — 12 mK; (a) Cb = Ct, (b) Cb = 1.25Ct- The inset is a detail of the structure indicated by *. 
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This figure "Figllanl.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0207437vl 



This figure "Fig21anl.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/0207437vl 



This figure "Fig31anl.GIF" is available in "GIF" format from: 
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